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Battery  alignment  is  one  o-f  the  critical  -factors  in 
the  control  o-f  the  weapons  systems  on  board  naval  warships. 
The  spatial  relations  o-f  the  guns,  directors,  and  launchers 
are  determined  by  re-ference  to  lines  and  planes  -fixed  with 
relation  to  the  hull  and  the  determination  o-f  these 
relations  constitutes  the  battery  alignment.  This  alignment 
only  compensates  -for  the  static  relations  between  the 
system  elements,  and  so,  when  these  planes  and  lines 
distort  in  a  seaway,  the  alignment  is  no  longer  correct. 

Normally,  the  re-ference  lines  used  -for  the  alignment 
are  the  ship  centerline  for  train  and  the  plane  of  the 
primary  element,  usually  the  director,  for  elevation.  Uhen 
the  alignment  is  done,  the  ship  is  alongside  a  pier  in  a 
known  condition  of  loading  and  in  calm  seas.  At  that  time 
the  relations,  vertical,  horizontal,  and  transverse, 
between  the  directors  and  the  weapons  are  measured  and 
recorded,  and  the  fire  control  system  is  adjusted  so  it 
compensates  for  static  parallax  when  making  the  fire 
control  solution.  Once  the  ship  goes  to  sea,  these 
relations  change  as  fuel  is  burned,  stores  used  and  the 
ship  moves  in  waves.  The  sea,  in  addition  to  causing  rigid 
body  motions  of  the  ship  hull  in  heave,  pitch,  surge,  sway, 
roll  and  yaw,  also  causes  it  to  flex  in  the  vertical  and 
horizontal  planes  and  to  twist  about  its  longitudinal  axis. 


This  -f1«xing  distorts  the  re-feramce  lines  established  in 
the  battery  alignment  and  introduces  errors  in  the 
targeting  solution  arrived  at  by  the  -fire  control  system. 
These  errors  are  due  to  dynamic  parallax,  resulting  -from 
relative  target  location  due  to  ship  attitude  and  -flexure. 

Studies  by  the  Naval  Sea  Systems  Command,  under 
contract  to  Rockwell  International  Ill,  have  shown  that  the 
-flexure  o-f  the  ship  hull  o-f  monohull  vessels  appears  to  be 
si  gn  i -f  i  cant ,  relative  to  rigid  body  motion.  It  seems  that  a 
method  to  measure  and  compensate  for  these  errors  is 
necessary  to  improve  the  fire  control  solution  and  so  to 
increase  the  fighting  effectiveness  of  the  ship. 

Measurement  of  the  flexure  could  be  done  directly,  for 
example  by  using  a  system  of  sighting  instruments,  and  the 
results  fed  to  the  fire  control  system  as  they  happen. 
Suitable  instruments  might  be  Lasers,  light  beams  or  strain 
gauges.  All  of  these  require  the  primary  equipment,  along 
with  its  support  systems,  which  are  subject  to  failure  and 
in  addition  the  system  must  be  aligned  and  kept  in 
alignment.  This  results  in  an  increased  maintenance  load  on 
ship's  personal,  and  the  reliability  problems  inherent  in 
mechanical  devices.  In  addition,  this  method  does  not 
provide  an  easy  capability  to  predict  what  the  flexure 
might  be  at  some  future  time,  thus  allowing  correction  of 
the  fire  control  solution  at  the  actual  instant  of  fire. 

An  alternative  method  involves  use  of  a  filter  to 


extract  de-flection  in-formation  -from  parameters  already 
measured  on  board  the  ship.  The  ship's  gyroscope  provides 
ship  heading  and  speed  in-formation  and  in  addition  it  can 
provide  heave,  pitch,  roll,  sway,  and  yaw  information. 
Information  regarding  wind  speed  and  direction  is  also 
available.  Additional  information,  if  needed,  might  be 
found  from  accelerometers  placed  at  the  bow,  the  stern  or 
other  critical  locations.  The  advantage  of  using  solid 
state  devices,  instead  of  directly  measuring  the 
deflections,  is  increased  reliability  and  descreased 
maintainace.  In  addition,  the  filter  formulation  lends 
itself  to  prediction  of  future  distortions,  which  can  be 
fed  to  the  fire  control  system. 

The  objective  of  this  project  is  the  formulation  of  a 
mathematical  model  of  the  ship/sea  system  that  can  be  used 
with  a  filter,  the  formulation  of  a  filter  to  estimate  the 
deflections,  and  a  predictor  to  estimate  future 
distortions.  The  equations  discribing  the  ship  can  be 
written  as  a  set  of  first  order  differential  equations, 
which  lend  themselves  to  a  state-space  model  formulation. 
It  is  then  natural  to  use  a  Kalman  filter  for  the 
estimator.  This  method  has  already  been  used  to  estimate 
the  rigid-body  motions  of  a  ship  for  input  to  an  automatic 
landing  system  [2,3,43. 
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CHAPTER  18  FORMULATIOj 


1.1  Literature  Ovepview 

Uave  induced  ship  vibration  is  the  subject  o-f  current 
interest  to  many  investigators.  The  beginnings  o-f  the 
majority  o-f  the  recent  work  dates  -from  the  late  1960's, 
when  the  per-formance  o-f  large  bulk  carriers  was 
investigated  by  Aertssen  tSl,  Aertssen  and  De  Lembre  [6], 
Bell  and  Taylor  171,  and  Goodman  [81.  These  papers  were 
mainly  concerned  with  the  statistical  nature  o-f  the 
stresses  and  bending  moments  and  their  prediction, 
Aertssen,  and  Bell  and  Taylor  doing  -full  scale  mesurements 
and  Goodman  taking  a  theoretical  approach. 

The  general  trend  in  the  theoretical  approaches  has 
been  to  use  a  modal  decomposition,  as  done  -for  example  by 
Bishop  and  Price,  et.  at.  19,10,11,12,133,  Chen  C14], 
Jensen  and  Pedersen  1153,  and  Kagawa  and  FuJ i ta  1163.  This 
has  an  advantage  when  working  with  regular  waves  and  for 
separating  out  the  contributions  of  the  various  modes. 

The  problem  of  the  exciting  hydrodynamic  forces  has 
been  addressed  by  Salvesen,  Tuck,  and  Fait  insen  [173  for 
the  rigid  body  motions,  using  strip  theory,  providing 
linear  expressions  for  the  distribution  of  the  wave  forces 
over  the  hull.  SkJ0rda1  and  Faltinsen  [183  have  developed  a 
linear  theory  for  the  hydrodynamics  of  a  springing  hull. 
Jensen  and  Pedersen  [153  provide  a  non-linear  formulation. 


B*tts,  Bishop  and  Price  tl93  tie  the  strip  theory-  approach 
together  with  structural  representations. 

The  use  o-f  state  space  techniques  -for  the  prediction 
of  hull  vibration  does  not  seem  to  have  been  used  by  any 
researchers  as  no  references  were  found  by  the  author 
di scribing  this  approach. 

1.2  Ship  Model 

There  is  no  established  procedure  on  how  to  best  model 
the  ship  and  the  sea  in  order  to  predict  flexural 
distortions.  One  major  consideration  is  the  size  of  the 
models  which  can  make  the  use  of  state-space  techniques 
very  cumbersome.  For  this  reason,  for  example,  it  is  best 
not  to  use  exact  structure  and  hydrodynamic  models,  such  as 
finite  element  techniques,  and  instead  to  try  to  reach  a 
compromise  between  model  size  and  modeling  accuracy.  A 
question  arises  when  trying  to  combine  a  structural  model, 
such  as  modal  representation,  with  a  random  sea.  The  sea, 
because  of  dispersion,  has  a  wavelength  which  is  a  function 
of  wave  frequency.  Even  though  the  random  sea  may  be  easily 
represented  by  Fourier  techniques  at  a  specific  point,  its 
representation  at  some  other  point  in  space  would  be  the 
result  of  the  propogation  of  an  infinite  number  of 

frequency  components.  How  to  include  these  phase  shift 
effects  in  a  modal  analysis  is  not  clear.  In  addition,  the 
computational  load  may  be  large  due  to  the  number  of  modes 
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required  to  eccuretely  represent  the  quasi-static  hull 
distortion,  given  that  the  natural  -frequencies  o-f  the  hull 
are  much  higher  than  the  wave  -frequencies  containing  the 
major  part  o-f  the  wave  energy. 

The  approach  used  here  is  to  approximate  the  hull  as  a 
quasi  static  non-uni -form  beam.  Using  simple  beam  theory,  one 
starts  by  considering  a  short  section,  dx,  o-f  the  hull  (see 
-figure  1.1).  The  governing  equation  o-f  the  beam  vibrating 
in  a  plane  o-f  symmetry  is  -found  by  summing  the  vertical 
forces  acting  at  point  "a* 

m<x)  dx  d/dt<dy/dt>  =  f<x,t)  dx  -  V<x,t>  *  V<x ,  t)■^dU/dx  dx 

=  f<x,t>  dx  +  (fJ/dx  dx  ( 1 . 1  > 


or 


m  d/dt<dy/dt)  *  f  V' 


( 1  .2) 


Similary,  summing  the  moments  about  'a*  one  obtains 

M'<x,t)  *  -U<x,t)  <1 .3) 


Then  by  using  (1.22)  in  (1.21)  one  finds 


m  d/dt<dy/dt)  **  f  -  M" 


<1  .4) 


From  simple,  linear  beam  theory,  one  has  the  relation 


M  =  Ely" 


I 


<1.5) 


Using  <1.24>  in  <1.23>  one  obtains 

m  d/dt<dy/dt)  =  -f  -  <EIy">"  <1.4)  , 

which  is  the  governing  equation  -for  vertical  vibrations  o-f 
a  beam.  ^ 

If  damping  is  allowed,  then  one  has 

I 

m  d/dt<dy/dt)  +  b  dy/dt  ♦  <EIy">"  =  f  <1.7> 


Once  one  has  the  governing  equation,  the  next  step  is  ) 

to  express  it  in  a  state-space  -form.  To  do  this,  let  xl  =  y 
and  x2  ®  dy/dt.  Then  with  this  coordinate  transformation, 

<1.26>  becomes  » 


dxl/dt  =  x2 

dx2/dt  =  -<EIxl")"  /  m-  bx2/m  +  f 


(1  .8) 


I 


The  next  step  is  to  represent  the  space  derivative. 
Using  a  finite  difference  approach,  one  lets  h  =  d;c ,  the 


I 
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I 


distance  between  points  on  the  hull.  Then  i -f  there  are  N 


points, 

L  *  <N-l)h  <1 .9) 

where  L  is  the  length  o-f  the  ship.  Expanding  the  space 
derivative,  one  has 

<EI  =  £1"  xi"  ♦  2  El'  xl'"  +  El  xi""  <1.18) 

Now  each  part  can  be  represented  seperately  by  -finite 
di -ff erences,  with  El,  El',  and  El"  as  input  data. 

At  the  ends  o-f  the  ship,  the  conditions  of  zero  shear 
and  zero  bending  moment  lead  to 


El  xl"  »  e 
<EI  xl")'  =  0 

■  El'  xl"  -►  El  xl'" 


<1  .11) 


Equations  (1.29)  with  boundary  conditions  (1.291)  can 
be  represented  to  the  order  desired  and  will  result  in  the 
formation  of  a  set  of  equations,  written  in  matrix  form  as 


3- 


where 


dXl/dt  =  A1*X1  ♦  B)«E 
Y  »  C1*X1 


(1 .12) 


X'  =  (xli ,xi2,x21 ,x22, . . . ,xNl ,xN2> 
xij  =  X  sub  J  at  point  i 

Y  s  vector  o-f  responses  at  the  N  points. 

A  =  a  2N  X  2N  matrix  containing  the  -finite 
di-f-ference  representation. 

E'  =  <fl,f2 . +N) 

B  »  distributes  the  ■fj'^s  to  the  proper  state 
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Oc*&n  waves  are  a  -function  o-f  both  the  wind  speed  and 
duration,  among  other  parameters.  One  o-f  the  commonly  used 
wave  spectra  that  attempt  to  account  -for  these  two  -factors 
is  the  Bretschne i der  Spectrum,  in  which  the  intensity  is 
speci-fied  by  the  signi-ficant  wave  height,  H,  and  the 
duration  by  the  modal  wave  -frequency,  Uhn.  This  spectrum  has 
the  -following  -form. 

S<w)  =  1.25  H««2  Uim»»4  exp<-l  .25<UJm/w)»*4>  <1.13> 

4  ws*5 

The  movement  of  a  ship  through  the  water  causes  it  to 
experience  excitation  not  at  te  wave  frequency,  w,  but  at 
the  encounter  frequency,  we,  where 

we  «  w  -  k’U  cosCbeta>  (1.14> 

The  energy  of  the  waves  remains  the  same,  however,  so 


S<w>dw 


S<we)dw« 


(1 .15) 


S(w>  dw  dw* 


th^rt'for* 


S<w*>  *  S<w> 
dw 
dw 


wh«p«  the  modulus  o-f  dwe/dM  is  used  because  S(we>  is  single 
sided. 

Tr i anta-fyll ou  and  Athans  C23  have  -found  a  rational 
approximation  to  (1.38)  in  the  -following  form 

S<we>  *  1  .25  H**2  6<a1pha> _ we«-»4  <1.17) 

4Wm  <  l-Kwe/wo)«»4)*#3 


(1 .18> 


alpha  s  u  ukn  cosCbeta> 

9 

B<alpha>  »  1.9339  <1.19> 

i*2  alpha 

wo(a1pha>  =  UJm<  l*a1  pha)  <1.20) 

0.8409 


The  corresponding  causal  system  is 


H<s)  *  So  _ ses2 _ 

< l+2J<s/wo)+<s/wo>**2>**3 


(1.21) 


where 


So  *  1 .25  H«»2  BCalpha)  (1.22) 

4  Mn 

J  «  0.707 

By  writing  (1.38)  in  standard  state  space  -form  <see 
chap  2)  one  has 


dXs/dt  =  As*Xs  *  Bs-ep 
ys 


Cs«Xs 


<1.23) 


where  p  is  wh 

i  te 

noise  o-f 

i  ntensi ty 

q  and 

e 

1 

0 

0 

0 

0 

:-wo*«2  - 

2Jwo 

0 

wo««2 

0 

0 

As  «  :  0 

0 

0 

1 

0 

0 

0 

0 

-wo**2 

-2Jwo 

0 

wo««2 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-wo**2 

-2Jwo 

Bs'«  [  0 

0 

0 

0 

0 

pi «wo«*3] 

Cs  *  C  So 

0 

0 

0 

0 

0  1 

< 1.24a) 

(1 .24b> 
( 1 . 24c ) 


Since  white  noise  hss  a  constant  spectrum  for  all 
f requenc i es 


Sp<w)  -  q 
P« 


<1 .25) 


Sp<w>  «  H(w>  *»2  S<w> 


wo**3  So  _ we»»4  wo»»4 

(wo*«4  +  we«s4>»«3 


(1 .24> 


one  may  idenitify  the  intensity  as 


q  «  p i  wo««3  So 


<1 .27) 
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1 .4  S9di  Model 


The  representation  o-f  the  sea  must  take  account  o-f  the 
dispersion  o-f  the  waves  in  space  and  still  allow  for  their 
random  nature.  If  the  relation  between  driving  point  1,  pi, 
and  distant  point  2,  p2,  is  dependent  on  a  fixed  wave 
number,  then  the  sea  at  p2  will  not  be  properly  related  to 
pi.  To  find  a  relation  that  will  give  the  correct  relation 
between  pi  and  p2,  one  starts  with  Lap1ace‘'s  equation  and 
the  free  surface  condition: 

Qxx  *  Qzz  »  0  in  the  fluid 

(1.29) 

Qtt  g  Qz  a  9  on  the  free  surface 


Then  by  assuming  that  the  potential  can  be  written  as 


Q(x,z,t)  =  U<x,t)  exp(kz>  (1.30) 


and  substituting  (1.41)  into  the  free  surface  condition, 
one  finds 

Qtt  gQz  =  (Utt  +  gkU)  exp(kz>  (1.31) 


If  the  deep  water  dispersion  relation,  k  =  w**2/g,  is  used, 
( 1 .42)  reduces  to 


Utt  +  w**2  U  *  0 


Utt  -  Utt  *  8 


shoMing  that  the  -free  sur-face  condition  is  satisfied  as 
expected.  Then  substituting  (1.41>  into  the  Laplace 
equation,  one  obtains 


g**2  Uxx  ♦  Utttt  =  0 


<1.32) 


It  is  intresting  to  note  that  <1.43>  may  be  obtained 
from  the  wave  equation  by  using  the  expression  for  wave 
phase  velocity  in  deep  water. 


c  =  w  /  k 


Wtt  -  c«*2  Wxx  =  0 


from  which 


UItt  -  w**2/<w**2/g)**2  Wxx  =  0 


w**2  Wtt  -  g**2  Wxx  »  0 
and  since  w»*2  W  =  -d/dt<dW/dt) 

Wtttt  ♦  g**2  Wxx  =  0 


Equation  <1.43)  can  now  be  represented  by  finite 
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di -f-f*rence  equations  in  space  and  as  a  set  o-f  -first  order 
equations  in  time  as 

dXv/dt  *  Av«Xw  *  Bw-»ys 

<1 .33) 

e  *  Cv*Xv  *  Dwsys 

where 

Xv'  a*  <ul  1  ,ul2,ul3|Ul4,u21 , . . .  ,uN4> 

uji  =  Uj  at  the  j-th  point 

uj2  *  dUj/dt 

uj3  =  d/dt<dUj/dt) 

uj4  =  d/dt<d/dt<dU4/dt>> 

E  »  vector  o-f  wave  elevations  at  the  N  points 


1.5  Force  Model 


The  hydrodynamic  exciting  forces  on  a  ship  hull  by  sea 
waves  have  been  studied  for  many  years.  The  approach  that 
seems  the  most  fruitful  for  this  project  is  the  strip 
theory.  Of  the  several  forms  available  [1^,193,  the  method 
of  Gerritsma  and  Beukelman  Ci93  will  be  used  here. 

The  hydrodynamic  exciting  force  on  the  ship  may  be 
di scribed  as  being  the  result  of  the  difference  between  the 
deflection  of  the  ship  hull  and  the  wave  profile,  z~y-w. 
Then  Gerritsma  and  Beukelman  have  proposed  that  the  force, 
F<x,t),  acting  on  the  ship  may  be  written  as 

F<x,t)  =  -  D  rMa<x)  Dz 1  ♦  N<x>  *  Dz  ♦  rgB<x>z  <1.34) 
Tt  L  DtJ  Ot 

where 

D  =  d  -  d  dx«d  -Ud 
dt  dx  dt  dt  dx 

Then,  using  z=y-w,  and  writing  the  force  as 


F<x  ,  t)  =  - 

•H<x 

,t)  +  Z<x,t> 

<1 .35) 

where 

X 

X 

=  lia(x)  D^y 

+ 

N<x)-U  dMa 

Dy 

+  rgB<x)y  < 1 .36) 

Dt* 

dx 

Dt 

Z<x,t> 

=  Ma<x)  D*w 

+ 

N<x>-U  ctfla 

Dw 

+  rgB<x)w  <1 .37) 

Ot*- 

dx 

5t 

the  force 

can  be  separated  into  the  part 

due  to  the  hu 1 1 

motions  and  the  part  due  to  the  waves.  The  part  due  to  the 
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hull  vibrations  may  be  moved  to  the  le-ft  hand  side  o-f  the 


equation  and  included  in  the  structural  model.  The  wave 
exciting  forces  then  may  be  written  in  a  form  that  may  be 
used  in  the  state-space  formulation.  Then  the  C-matrix  and 
D-matrix  of  the  sea  model  are  modified  with  equation  <1.53> 
to  output  the  force  instead  of  the  elevation.  The  output 
matrix  of  the  spectral  model  is  modified  to  output  the 
first  and  second  time  derivatives  of  the  wave  elevation  in 
addition  to  the  elevation  itself. 

1.6  Total  model 

The  models  for  the  sea  spectrum,  the  sea  surface,  and 
the  ship  can  be  combined  into  one  model  by  identifying  the 
output  of  one  with  the  input  of  another  in  cascade  fashion. 


where 


dX/dt  *  A*X  +  B«p 
Y  »  C*X 


(1.38) 


A 

B' 

C 

X" 


P 


:  Ai 

Bl*Cv 

Bl*Dv«Cs 

:  0 

Av 

Bv*Cs 

:  e 

0 

As 

C  0 

0 

Bs' 

[  Cs 

0 

0 

<  XI  ' 

Xv' 

Xs' 

N<0,1) 
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CWiPTER  2;  STATE  SPACE  TECHNIQUES 


2.1  State  Space  Models 

In  order  to  apply  Kalman  Filtering  techniques,  one 
must  -formulate  the  problem  in  the  -form  o-f  a  state-space 
model.  This  is  a  set  o-f  -first  order  di -f-ferent  i  a1  equations 
that  di scribe  the  dynamics  of  the  system,  its  inputs,  and 
i ts  outputs. 

The  state  of  a  system  may  be  defined  at  any  time  as  a 
minimum  <not  neccessarily  unique)  set  of  parameters  which 
must  be  specified  to  di scribe  that  system.  A  state  vector 
is  a  vector  composed  of  these  parameters  (state  variables). 
The  number  of  state  variables  required  is  equal  to  the 
order  of  the  system,  n. 

For  example,  a  simple  mass-spring-damper  system  is 
discribed  by  a  second  order  differential  equation. 

m*d/dt<dz<t)/dt)  -*■  b*dz<t)/dt  *  l<*z<t)  =  f<t)  (2.1) 

Then  one  may  define  the  state  variables  as  the 
position,  z,  and  the  velocity,  dz/dt,  of  the  mass.  Next  let 
xl  s  z  and  x2  »  dz/dt.  Uli  th  this  change  of  variables,  (2.1) 
may  be  written  as 

dxl/dt  =  x2 


dx2/dt  *  -k/m*xl 


b/m*x2  f/m 


(2.2) 


or  in  matrix  notation 


where 


dX/dt  *  A*X  +  B»F 


A  = 

B'  = 
X'  = 


:  e  1  : 

;  -k/m  -b/m  ! 
C  e  1  ] 

(  xl  x2  ) 


<2. 3) 


Alternatively,  because  the  state  variables  are  not 
unique,  one  may  de-fine  the  state  variables  as  the  -force  in 
the  spring,  K,  and  the  momuntum  o-f  the  mass,  P.  A 
trans-formation  may  be  found  from  one  to  state  vector  to  the 
other  by  using 


K  «  k*z 
P  =  m*dz/dt 


(2.4) 


One  may  also  define  an  input  vector  U,  containing  the 
r  inputs  to  the  system  and  an  output  vector,  Y,  containing 
the  m  outputs.  The  system  equations  are  then  a  set  of  n 
first  order  differential  equations 


dxi/dt  »  fl<X,U,t) 


<2. 5) 


dxn/dt  *  fn<X,U,t) 


dX/dt  =  F<X,U,t) 


<2.^) 


The  output  equations  relate  the  output  vector,  Y,  to 
the  state  vector  and  the  input  vector 


yl  =  gl<X,U,t> 


ym  =  9m<X,U,t) 


Y  =  6<X,U,t) 


(2.7> 
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If  the  system  is  linear  or  linearized,  one  may  write 
the  functions  fi  and  gi  as  linear  combinations  of  the 
states  and  the  inputs. 


<2.8> 


<2.9> 


dX/dt  -  A*X  ♦  B*U 
Y  =  C*X  +  D*U 


(2.10) 


where  A  i s  an  n  x  n  matrix  di scribing  the  system 
dynamics,  B  is  n  x  r  and  di scribes  the  distribution  of 
inputs,  C  i s  an  m  X  n  matrix  discribing  the  distribution  of 
states  to  the  outputs,  and  0  is  a  m  x  r  matrix  discribing 
the  distribution  of  inputs  to  outputs.  The  form  di scribed 
by  (2.10)  is  known  as  Standard  State-space  Form. 


2.2  Kalman  Fi 1  ter 

In  this  problem  the  object  is  to  be  able  to  estimate 
the  flexural  response  of  the  ship  hull  at  time  t  and  from 
this  to  predict  what  the  response  will  be  at  time  t+dt.  To 


do  this,  one  common  method  is  to  use  the  Kalman  Filter. 
This  is  a  method  which  derives  an  optimal  (in  the  least 
squares  sense)  estimator  -for  the  state  o-f  a  system 
represented  by  a  state-space  model . 

Assume  one  has  a  mathematical  model  o-f  the  system 
under  consideration, 

dx/dt  *  A*x  *  B*w 

(2.11) 

P  =  C*x 

where  w  is  white  noise  of  intensity  Q,  the  expected  value 
of  X  at  ta8,  ECx(8)],  is  x8  =  8,  and  the  covariance  of  x  at 
t*8 ,  ECx(8)x^ (8) ] ,  is  X8 .  Assume  further  that  one  has  1 
measurements,  z,  which  are  related  to  x  by 

z  =  M*x  V  (2.12) 

where  v  is  measurement  noise,  modeled  as  white  noise  of 
intensity  R,  uncorrelated  with  the  process  noise,  w. 

Also  assume  that  the  system  is  observable,  meaning 
that  the  initial  condition  of  any  state  may  be  found  by 
measuring  the  output  for  a  finite  period  of  time. 

Then  one  may  propose  an  observer  of  the  form 

dy/dt  =  A*y  -►  K*<  z  -  M»y  )  (2.13) 

y<8)  =  E[  x<8>  3  =  8 

K  is  known  as  the  Kalman  gain  matrix. 


Ont  may  than  de-fine  the  error  in  the  state  estimate  as 


•  -  y  -  X  <2.14) 

and  the  error  covariance  as 

P  »  EC  ee'  3  <2. 15) 

It  may  be  shown  [231  that  the  equation  which  controls 
propagation  o-f  the  error  covariance  is 


dP/dt  *  <A-KM)P  ♦  P<A-KM)'  -I-  KRK'  BQB'  <2.16) 

P<0)  *  X<0)  =X0 

Examining  the  covariance  matrix,  P,  finds  that  the 
diagonal  is  made  up  of  the  variances  of  each  of  the  states, 


P  =  EC  ee' 


•  2  ] 


en 


el««2  el«e2  . 
e2«el  e2««2  . 

e  s  • 


el  e2  ...  en 


en**2  1 


<2.17) 


So  if  one  takes  he  Trace  of  P 


Tr<P> 


n 


e  i  **2 


<2.18> 


one  has  the  sum  o-f  the  squares  o-f  the  errors  in  each  o-f  the 
states.  Athans  C23]  shows,  using  the  matrix  minimum 
prinicipal,  that  the  K  that  minimizes  the  Trace  o-f  P,  ie; 
the  sum  o-f  the  squares  o-f  the  error,  is 


K  =  PM'R  <2. IP) 

I-f,  as  in  the  present  case,  none  o-f  the  matrices  A,B,M,Q 
and  R  are  -functions  o-f  time,  then  dP/dt  -^0  in  the  steady 
state  solution,  and  one  has 

e  at  AP  4  PA^  4  BOB'  -  PM'R  MP  <2. 20) 

Summarizing,  then,  i-f  one  has  a  linear  system 
described  by 


dx/dt  ■  A-ex  4  B«w 
p  *  Cx 


(2.21) 


one  may  estimate,  with  minimum  square  error,  the  state  x  by 
using  y  which  is  described  by  the  equation 


dy/dt  «  A*y  4  K<z-My) 
K  =  PM'R 


<2.22) 

<2.23) 
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with  P  giv*n  by  (in  the  steady  state  case) 

0  »  AP  ♦  PA'  +  BQB'  -  PM'R  MP  (2.24) 

So-ftware  packages  exist,  -for  example  CTRLC  [26],  which  will 
solve  2.24  -for  the  steady  state  values  in  P  and  so  -for  the 
values  in  K. 

Once  the  values  o-f  the  state  at  time  t0  are  known, 
either  -from  direct  measurement  or  from  estimates,  the  r  -xt 
task  is  to  predict  the  the  states  at  some  future  time, 
tf»t0-»dt.  To  do  this,  one  may  solve  2.291  subject  to 
x(t0)«x0  as 


x<t)  *  exp(A< t-t0) >»x0  * 


expCA(t-s)>Bw(s>ds 

't0 


(2.25) 


Now,  suppose  that  at  time  t0,  one  wants  to  predict 
x<  t0-*-dt)  .  Then 


tf 


x<tf)  ■  expCACdt)  >*x0  -♦■ 


expCA(tf-s)>6w(s>ds 

t0 


<2. 26) 


but  t0<s<atf  is  some  future  time  and  w(s)  is  unknown 
because  w  is  random.  In  fact,  since  w  is  white  noise  and  so 
completely  random,  it  is  completely  unknown  and  the  est 
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estimate  o-f  it  is  its  average  value  or  expected  val-ue, 


E[  w(t>  ]  »  8  (2.27) 

There-fore,  the  best  estimate  o-f  x<t-f>  is 

xp<t-f>  * 

where  x  has  been 
is  obtained  -from 


exp<Adt>»x<te>  <2.28> 

used  in  the  initial  condition  because  it 
the  Kalman  -filter. 
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2.3  Approximations 

Th*  question  o-f  the  order  o-f  the  'finite  difference 
representation  will  be  addressed  next.  The  accurarcy  of  the 
distortion  results  are  dependent  on  that  order.  For  just 
one  point  of  the  ship  model,  say  J,  equation  1.27  can  be 
written  as  (with  J  understood  in  the  coefficients) 


dxjl/dt  =  xj2 

dxj2/dt  =  xjl"" 

m 

xjl" 

m 


-  El'  xjl"'  - 
m 

-  ^  xJ2  +  _f^ 

m  m 


(2.29) 


Since  the  highest  derivative  in  x  is  the  controlling  factor 
in  the  discretation  error,  the  case  of  a  uniform  ship,  with 
El'  and  El"  both  equal  to  rero  will  be  investigated.  The 
method  used  will  be  that  used  by  ^i  th  [233. 

At  point  j,  (2.30)  becomes 


dX/dt  =  A*X""  +  D*X  +  F 


(2.30) 


where  the  matrix  A  is 


and  D  is 

D  » 


0  0 

-El/m  0 

0  0 

0  -b/m 
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Now  write  the  -fourth  derivative  o-f  X  as 


«  X<J-»2)  ~4X<j-H)  ♦  6X<j>  -  4X<j~l)  ♦  X<j-2) 

h**4 

Using  (2.32)  in  (2.31),  one  -finds 

dX<j)/dt  =  A  X<J->-2)  -4AX<J-H)  ♦  6A  X(J> 

h**4 

-  4<i>X<J-l)  ♦  AX<j-2) 
h**4 

*  0  X<j) 

and  expanding  the  X<j-*-l)  by  Taylor  series 

X<j-H)  *  X  *  hex"  J__h**2*X"  *  l__h**3*X" '  *  ... 

2!  3! 

Substituting  (2.34)  into  (2.33),  Keeping  terms  to 
h««8  and  collecting  terms  leads  to 

dX<j>/dt  »  A»X""  *  D»X  +  F  * 

<  JL»A*X®  *h**2  )  *  <  _l_*A*)C®°^h«*4  ) 

6  80 


(2.31) 


(2.32) 


(2.33) 

order 


(2.34) 


where  the  term  in  the  -first  bracket  can  be  interpeted  as 
the  order  h**2  error  and  that  in  the  second  bracket  as  the 
order  h»»4  error. 

In  order  to  see  the  e-f-fect  on  a  particular  solution  to 
(2.30)  y  consider  the  solutions  o-f  the  natural  modes  o-f  a 
uni -form  beam  [223  .  The  space  dependent  part  is  composed  o-f 
sinusoids  and  hyperbolic  -functions  in  the  variable  kx. 

w<x>  =  sin<kx>  +  sinh<kx>  +  s»<  cosCkx)  *  cosh(kx)  )  <2.35) 

where  s  is  a  constant  dependent  on  the  mode.l-f  L  is  the 
length  then 

kL  <4*  )*Pi 

de-fines  the  approximate  values  o-f  k.  Taking  the  eight 
dirivative  o-f  <2. 35)  brings  out  a  k-e^S  term,  so 

w  ■  k««8«w 

and  <2.31)  can  be  written  as 


dX/dt  ■  A*X"" 


*  D*X  -»  F  E 


(2.36) 


wh«re 


E  =  A*X*°^ *h**4 

88 

=  A*k*«8*X*h**4 

88 


Then  since 


k*«4*X  =  X"" 

one  may  wr i te 

dX/dt  =  A*X''"  <  1  +  1  sh»»4»k«4  )  +  D*X  +  F 

88 

Then  the  relation  between  the  point  spacing  and  the 
may  be  found  from  the  second  term  in  the  brackets 

e  =  J _ *h«*4*k**4 

88 

If  the  maximum  desired  error  is  eo,  then 

h**4  <  88*  eo 


<2.37> 


error 


k»*4 


or 


h  <  2.99*L  eo 

<j+i)Pi 


<2.38) 


By  similar  means  the  relation  -for  order  h-squared  can 
be  -found  to  be 


h  <  2.45«L  eo  <2. 39) 

TJ+J)PT” 

Table  2.1  shows  h  vs.  eo  -for  both  orders  -for  the 
sprinQing  mode  f  and  table  2.2  shows  h  vs.  j  -for  a  -fixed 
eoa>6.12S.  Length  is  100  meters. 


Table  2.1 
h  vs.  eo  < J*1 ) 


h**2 


h**4 


eo 

h 

"Npts 

h 

Np  ts 

0.5 

36.76 

3 

53.35 

2 

0.25 

26.00 

4 

44.87 

3 

0.125 

18.36 

6 

37.73 

3 

0.0625 

13.00 

8 

31.72 

4 

Table  2.2 


h  vs.  j  <eo  =  8.125) 

h«»2  h«»4 

J  h  Npts  h  Npts 

1  18.38  6  37.73  3 

2  11.03  9  22.64  5 

3  7.88  13  16.17  7 

4  6.13  17  12.58  8 

The  most  obvious  -fact  to  be  gained  -from  Table  2.1  is 
that  -for  a  given  accuracy,  the  number  o-f  points  and  hence 
the  number  o-f  states  is  much  less  with  the  order  h')^'K4 
approximation.  Each  added  point  requires  six  additional 
states  in  the  total  model,  two  -for  the  ship  model  and  four 
for  the  sea  model.  However,  only  five  points  are  required 
to  approximate  the  fourth  direvative  to  order  h*-*2  but 
eight  points  are  required  for  order  h**4.  These  numbers  set 
the  minimum  number  of  points  required  for  any  accuracy. 

As  the  number  of  nodes  in  the  vibration  increases,  the 
number  of  points  for  a  given  accuracy  increases  rapidly,  as 
seen  in  Table  2.2.  Again,  the  order  h«'»4  approximation  is 
the  less  restrictive  of  the  two.  It  only  requires  8  points 
to  achieve  12.5>i  accuracy  with  order  h««4,  while  17  points 
are  required  for  the  order  h««2  approximation. 

This  means  that  while  it  is  possible  to  use  an  order 
h««2  approximation,  an  actual  i npl ementat i on  will  require 
both  a  higher  order  and  a  greater  number  of  points. 


CHAPTER  3;  NUMERICAL  RESULTS 

3.1  Stability  o-f  Models 

The  mathematical  models  o-f  the  sea  spectrum  and  the 
ship  structure  are  stable  as  modeled.  However  the  -force 
model  is  not.  The  -fourth  order  equation  di  scribing  the  sea 
sur-face  has  -four  solutions,  two  o-f  which  are  exponentials, 
and  it  is  these  solutions  that  can  cause  instability.  The 
equation  is 


Utttt  *  9**2  Uxx  =  0  <3.1) 

which  has  solutions  in  the  -form 
U<x,t)  ~  C  exp<iwt),  exp<-iwt>,  exp(wt>,  exp<“wt)  > 

<3. 2) 

The  exponential  roots,  in  the  analytical  problem,  are  taken 
care  o-f  by  the  initial  conditions.  UJhen  using  a  computer, 
however,  even  i -f  the  initial  conditions  are  set  equal  to 
zero  subsequent  round-o-f-f  errors  introduce  -fictitious 
non-zero  initial  conditions,  which  excite  growing 
exponent i a1 s. 

To  see  this,  consider  the  simple  example 

Ytt  -  a**2  Y  =  0  <3.3> 

Y<0)  a  0 

Y< i nf i n i ty)  *  0 

which  has  the  analytic  solution 


Y<t)  a  A  exp(at)  -»■  B  exp<-at) 


<3.4> 


By  approximating  this  by  -finite  di -f -f erences ,  one  obtains 


Ytt  ~  Y(i-H)  -2»Y<i)  *  Y<i~l)  =  a»*2  Y<  i  )  <3. 5) 

h**2 
or 

Y<i-H)  -<2  +  <ah)**2)  Y<  i )  ♦  Y<i-1>  »  8 
To  solve  this  -finite  di-f-ference  equation,  take 

Y<n)  =  r**n 

then  one  obtains 

r**2  -<2  -*-  <ah)**2)  r  -f  1  =  8 
and  then  -for  h  <  1/a 

r  ~  1  -►  ah  +  0<ah**2) 

and 

r  ~  1  -  ah  0<ah*«2> 

Then  taking  a  linear  combination  o-f  the  solutions 

Y<n)  =  A  <l  -f  ah)**n  -♦■  B  <  1  -  ah>*»n  <3. 6) 

it  can  be  seen  that  no  matter  how  small  A  is,  ii  it  is  not 
exactly  equal  to  zero,  the  part  o-f  the  solution  involving 
the  growing  exponenitial  will  eventually  dominate  the 
solution,  swamping  the  correct  solution.  I-f  the  objective 
were  only  to  obtain  the  solution  to  the  equation,  one  might 
try  integrating  backwards  in  time.  Here  though,  the 
objective  is  to  produce  the  solution  in  real  time,  so  the 
equation  must  be  integrated  forward  in  time. 

This  is  not  as  severe  a  problem  as  it  might  seem  at  first 
look.  One  of  the  properties  of  the  Kalman  Filter  is  the 
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stabil  i  zat  i  on  o-f  the  -filter  model.  The  gain  that  results 
will  result  in  a  closed  loop  observer  matrix, 

Acl  «  A  -  K*C 

that  has  all  its  eigenvalues  in  the  le-ft-hal-f  plane.  The 
only  place  the  problem  arises  is  in  simulation  o-f  the 
system.  By  de-fining  a  reconstruction  error,  e=x-xhat ,  and 
using  the  expressions  -for  x  and  xhat,  it  can  be  seen  that 
the  error  must  obey  the  di -f-f  erent  i  a1  equation 

de<t)  =  CA-KC]  e<t)  <3.7> 
dt 

This  means  that  the  reconstruction  error  goes  to  zero  as  t 
becomes  large  only  i-f  the  closed  loop  observer  m»trix  is 
asymptotically  stable. 

Simulation  o-f  the  system  is  necessary  to  test  the 
observer  as  the  model  is  varied.  The  method  used  is  to  sum 
a  number  o-f  sinusoids,  each  o-f  which  has  an  amplitude 
proportional  to  the  energy  in  the  sea  spectra  at  that  wave 
-frequency  and  a  phase  that  is  random  with  respect  to  the 
phases  o-f  the  other  sinusoids. 

Assume  that  we  are  given  a  sea  spectrum,  S<w) ,  and  we 
wish  to  simulate  the  sea  elevation  at  a  number  o-f  points. 
Since  the  instantanous  sea  elevation  at  a  point  x  and 
frequency  w  can  be  represented  as  the  square  root  of  the 
energy  at  that  frequency  (multiplied  by  2  since  the  sea 


spectra  is  only  de-fined  -from  8  to  infinity), 

V<x,w>  =  V2  S<x,w)’ 

then  the  elevation  at  x,  including  all  frequencies  should 
be 


oo 


(3.8) 


where  p  is  the  random  phase.  Then  including  the  time  factor 
one  has 


lOO 


V<x)  =  Re<  \12  S<w)  dw'  expf  i  <wt  *  p<w)  >)  (3.9) 

e 


For  a  computer  simulation,  equation  <3. 9)  must  be 
discretized  to  become 


V<x,t> 


Nw 


2  S<wJ>  cost  wj  t  p(wj)  )  (3.10) 


j=l 


where  here  p  is  a  random  number  uniformly  distributed 
between  0  and  two  pi. 


The  simulated  sea  surface  is  then  used  to  drive  the 


ship  model  to  -form  s  simulation  model,  the  output  o-f  which 
is  used  as  a  source  o-f  the  measurements  supplied  to  the 
observer  model.  This  model  was  run  and  the  results  stored 
in  a  computer  -file  -for  later  use.  Figure  3.1  shows  a 
representative  run  o-f  the  simulation. 

3.2  Estimation  and  Prediction  with  Correct  Model 

The  Kalman  filter  for  the  total  model  is  required  to 
construct  optimal  estimates  of  the  states  of  the  system 
that  can  be  used  to  predict  the  future  output  of  the 
system.  This  includes  those  states  that  cannot  be  measured 
directly,  such  as  those  of  the  sea.  To  do  this  it  must  have 
a  source  of  information  about  the  system. 

The  filter  gains  are  affected  not  only  by  the  model 
used  but  also  by  the  noise  in  the  measurements  themselves. 
This  noise  results  from  a  number  of  sources.  For  example, 
from  measuring  instrument  noise,  slamming  of  the  hull, 
local  vibration  effects  and  vibrations  excited  by  external 
sources.  The  filter  then  must  reject  the  higher  frequency 
noise  to  obtain  the  needed  information  about  the  states. 
The  measurements  are  found  from 

z  «  y  V,  V  "  N<0,R)  (3.11> 

where  R  is  the  measurement  noise  covariance  matrix,  defined 
as  the  expected  value  of  the  noise  process 
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Figure  3.1 

Simulated  Hull  Deflection  at 
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T  I  ne 


Then  what  is  needed  are  estimates  o-f  the  measurement  noises 
that  would  be  encountered  in  an  actual  ship  hull.  Using 
studies  o-f  ship  ui  brat  ion,  the  noise  covariance  for 
displacement  was  taken  as 

r(i|i>  31  3.8  E  -4  meters  squared 

at  amidships.  According  to  Lewis  127],  the  amplitude 
amidships  will  be  equal  to  about  8.25  the  value  at  the  end 
for  machinery  induced  vibration.  For  velocity  of 
displacement  a  covariance  of 

r<i,i>  3  6.8  E  -1  meters  squared  /  second  squared 

at  amidships  was  used.  All  off  diagonal  values  were  assumed 
to  be  equal  to  zero,  implying  that  the  measurement  points 
are  isolated  from  each  other. 

The  filter  was  calculated  using  the  models  for  the  sea 
spectra,  forces  and  ship  in  cascade.  This  allowed  an 
investigation  to  determine  the  required  measurements  to 
insure  the  observability  of  the  system.  In  addition,  the 
filter  pole  locations  were  determined.  For  the  Great  Lakes 
Ore  Carrier  Stewart  J.  Cort,  described  in  Appendix  A,  the 
poles  were  all  located  at  a  frequency  of  less  than  1.85 
radians  per  second.  The  pole  locations  are  shown  in  Table 
3.1  for  the  output  matrix 


:  1  e  e  e  e  8  0  e 
C»  81088888 
0  0  1  8  8  8  8  8 
8  8  8  1  8  8  8  8 

This  output  mstrix  provides  measurements  o-f  the 

displacement.  The  resulting  displacement  estimation  errors  * 

are  given  in  Table  3.2  and  Figures  3.2  and  3.3  show  the 

estimated  value  compared  with  the  simulated  value.  The 

noisy  signal  is  shown  in  -figure  3.4.  * 


0.03 


Figure  3.2 


Estimated  vs.  Simulated 
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Figure  3.2 
Noisy  Measurement 
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RMS  Estimation  errors 
<tneters) 


9.54E-2 
9.3<4E--2 
1 .54E-1 
9.50E-2 


Prediction  o-f  the  -future  values  o-f  the  states,  given 
the  current  measured  or  estimated  values  can  be  done  by  the  • 

equat i on 


dxp  »  A  *  xp  (3.14)  • 

dt 

xp<e>  »  xhat<e> 

I-f  A  is  stable,  the  predicted  value  will  go  to  zero  as  time  _• 

gets  large.  The  objective  then  is  to  determine  how  -far  in 
the  -future  the  prediction  is  acceptable.  If  dt  is  the 
prediction  time  step,  the  actual  prediction  becomes 


50- 


xp<t)  »  expC  A  *dt  >  *  xhat<t>  <3.15) 

which  involves  matrix-vector  multiplication  only. 

For  this  system,  a  number  o-f  prediction  models  were 
tried,  but  none  gave  usable  predictions.  The  -first  model 
tried  was  to  use  the  unstabilized  system,  but  as  expected, 
the  output  immediately  started  growing. 

Next,  the  system  model  was  stabilized  by  using  a  very 
high  measurement  noise  to  calculate  the  -feedback  gain.  For 
this  model,  the  predicted  de-flections  showed  an  initial 
change  in  the  same  direction  as  the  true  value  but  rapidly 
decreased  in  amplitude  to  about  one  tenth  o-f  its  amplitude. 
The  phase  o-f  the  prediction  was  very  close  to  that  o-f  the 
true  de-flection. 

Another  method  tried  was  to  stabilize  the  model  of  the 
sea  surface  only.  This  resulted  in  a  stable  total  model  but 
no  improvement  was  seen  in  the  prediction.  Here  the 
predicted  value  rapidly  diverged  from  the  true  value  after 
about  four  time  steps  and  the  phase  showed  no  apparent 
relation  to  that  of  the  true  phase. 

Finally,  a  prediction  was  tried  by  using  only  the 
model  of  the  ship.  This  was  even  less  successful  as  the 
prediction  immediately  diverged  from  the  true,  often 
increasing  when  the  true  value  was  decreasing  and 
V i ce-versa. 

Though  prediction  was  not  acheived  by  any  of  these 
models,  this  does  not  mean  that  prediction  can  not  be  done. 
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It  is  one  o-f  the  underlying  ideas  o-f  the  state  space 
technique  that  all  necessary  ih-format  i  on  about  a  system  is 
included  in  the  state  o-f  the  system  at  any  time.  Prediction 
models  that  bare  investigation  are:  (1>  stabilize  the  sea 
sur-face  model  by  treating  it  as  a  regulator  problem  and 
selectively  choosing  the  elements  in  the  control  weighting 
matrix  to  limit  the  wave  amplitude,  (2>  try  the  same  method 
on  the  total  model,  <3>  a  useable,  i-f  short  prediction  may 
result  -from  the  second  method  tried  i-f  the  time  step  is 
reduced,  <4)  use  the  -first  method  tried  but  scale  the 
predicted  output. 
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CHAPTER  4;  DISCUSS I ON 


The  problem  addressed  in  this  thesis  is  the  estimation 
o-f  ship  hull  vibrations  -from  a  minimum  number  o-f  measured 
parameters.  These  estimations  are  then  used  to  predict  the 
future  deflections  of  the  ship  hull. 

The  first  step  was  to  develop  a  mathematical  model  of 
the  ship-sea  system  in  state-space  form.  This  involved 
three  separate  models.  The  first,  which  models  the  sea 
power  spectrum,  was  taken  from  the  literature.  This  model 
is  based  on  the  ITTC  ocean  spectrum  that  uses  as  parameters 
the  modal  wave  frequency  and  the  significant  wave  height. 
In  this  spectra,  the  modal  frequency  can  be  taken  as 
4.47/Tave,  where  Tave  is  the  average  period  measured  at  the 
calm  water  level.  The  significant  wave  height  may  be  taken 
from  table  4.1  i f  no  measurements  are  available. 


wind  speed  (knots) 
20 
30 
40 
50 
60 


sig.  wave  height  (m> 

- 571 - 

5.1 

8.5 

11.0 

14.6 


The  model  to  include  the  dispersion  relation  for  sea 
waves  was  developed  from  a  velocity  potential 
representation  of  the  sea,  resulting  in  an  unstaolc  model. 
The  instability  was  due  to  the  increasing  exponential 
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solution  to  the  resulting  -fourth  order  in  time  di  f-ferent  i  al 
equation. 

The  third  mode]  models  the  ship  hull  itself.  This  was 
derived  by  assuming  that  the  ship  could  be  represented  as  a 
simple  beam.  Then  a  composite,  continuous  time,  discrete 
space,  finite  difference  representation  of  the  simple  beam 
equation  was  used  for  the  final  model.  The  model  was 
developed  using  a  second  order  finite  difference 
approximation  but  for  an  actual  application,  a  fourth  order 
approximation  should  improve  the  accuracy  to  a  significant 
degree . 

Gerritsma  and  Beuke1man''s  relation  relating  the  wave 
height  to  force  on  the  ship  was  used  to  couple  the  sea 
model  with  the  hull  model.  Since  in  a  state-space 
formulation  the  time  derivatives  of  the  wave  height  and  the 
deflection  are  available  in  the  state  vector,  implementing 
the  force  model  simply  required  additional  terms  in  the 
output  matrices  of  the  sea  spectra  and  sea  surface  models, 
and  added  terms  in  the  hull  system  matrix. 

Using  these  models,  CTRLC  C26],  a  software  package  for 
control  system  design,  was  used  on  a  Digital  k^AX-1 1/782  at 
the  Joint  Computer  Facility  to  do  the  simulations,  develop 
the  filter  gains  and  evaluate  the  sensitivity  to  errors  in 
the  model.  The  estimation  of  the  deflections  of  the  hull 
could  be  done  well  with  the  state  space  model  developed. 
The  prediction  of  the  deflection  to  some  time  in  the  future 
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was  not  succsss-ful  using  th*  models  developed  here,  but 
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APPENDIX 


PREDICTION  RESULTS 

The  several  methods  o-f  prediction  attempted  gave  less 
than  sat i s-factor/  results  when  compared  to  the  simulated 
signal.  Figure  A.l  shows  the  prediction  that  resulted  -from 
using  the  system  model  stabilized  by  using  a  very  large 
measurement  noise.  This  resulted  in  a  predicted  value  that 
was  much  larger  than  the  simulated  value.  The  same  general 
pattern  resulted  when  the  measurement  noise  was  varied  -from 
a  low  equal  to  the  expected  measurement  noise  to  as  high  as 
could  be  computed.  Nor  did  changing  the  starting  point  o-f 
the  prediction  with  relation  to  the  phase  o-f  the  simulated 
signal  improve  the  prediction. 

The  next  -figures  show  the  prediction  resulting  -from  a 
model  constructed  by  stabilizing  the  sea-sur-face  model 
be-fore  -forming  the  total  model.  In  Figure  A. 2,  the 
predicted  signal  is  very  similar  to  the  prediction  in 
Figure  A.l.  I-f  Just  the  -first  -few  time  steps  are  plotted, 
it  can  be  seen  that  the  predicted  signal  does  not  track  the 
simulated  signal  even  -for  a  short  period  o-f  time,  Figure 
A. 3. 

Finally,  Figure  A. 4  shows  the  results  o-f  using  only 
the  ship  model  to  do  the  prediction.  As  can  be  seen,  this 
prediction  has  no  apparent  relation  to  the  simulated 


